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O ■ Abstract. 

We present a review on the determination of the primordial helium abun- 
dance, Y p , based on the study of hydrogen and helium recombination lines in 
£^ ' extragalactic H n regions. We also discuss the observational determinations of 

d [ the increase of helium to the increase of oxygen by mass AY/ AO, and compare 

them with predictions based on models of galactic chemical evolution. 
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> ■ 1. Why Y p ? 

m : 

The determination of Y p is important for at least the following reasons: a) 
It is one of the pillars of Big Bang cosmology and an accurate determination 
of Y p permits to test the standard Big Bang nucleosynthesis (SBBN), b) the 
t — , models of stellar evolution require an accurate initial Y value; this is given by 

Y p plus the additional Y produced by galactic chemical evolution, which can be 
estimated based on the AY /AO ratio, c) the combination of Y p and AY/ AO 
Q-f is needed to test models of galactic chemical evolution, d) to test solar models 

it is necessary to know the initial solar abundances, which are different to the 
photospheric ones due to diffusive settling, this effect reduces the helium and 
heavy element abundances in the solar photosphere relative to that of hydrogen, 
the initial solar abundances can be provided by models of galactic chemical 
evolution, e) the determination of the Y value in metal poor H n regions requires 
a deep knowledge of their physical conditions, in particular the Y determination 
depends to a significant degree on their density and temperature distribution, 
therefore accurate Y determinations combined with the assumption of SBBN 
provide a constraint on the density and temperature structure of H n regions. 
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2. Recent Determinations of Y„ 



Previo us r eviews on Y p det erminations have been presented bv lPeimbert et al 



2003) a ndlLuridianal (120031) . Re cent determinations of Y p are those bvllzotov et al 
ll99Sh. Ilzotov fc Thuanl (12004T). Ilzotov et al.1 (1200611 . iLuridiana et al.1 (120031) " 



Olive fc Skillmanl (|2004l ). lFukugita fc Kawasakil 020061 ) ■ and lPeimbert et al.1 (|2007h 
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A crit ical discussion of these determinations has been presented bv lPeimbert et aT 
(120071 ). Most of the differences among these determinations are due to system- 
atic effects. 



3. Error Budget 

The error budget of the Y p determination bv lPeimbert et aT ( 20071 ) is presented 
in Table 1. In this table the sources of error are listed in order of importance. 
Th e error budgets of Y p determinations by other groups are different to that 
by iPeimbert et aD ((2007) for many reasons, each error budget depends on the 
sample of H n regions used and on the treatment given to the different sources 
of error. 



Table 1. Error budget in the Y p determination 



Problem Estimated error 



Collisional Excitation of the H I Lines 


±0 


0015 


Temperature Structure 


±0 


0010 


O {AY/ AO) Correction 


±0 


0010 


Recombination Coefficients of the He I Lines 


±0 


0010 


Density Structure 


±0 


0007 


Underlying Absorption in the He I Lines 


±0 


0007 


Recombination Coefficients of the H I Lines 


±0 


0005 


Underlying Absorption in the H I Lines 


±0 


0005 


Ionization Structure 


±0 


0005 


Collisional Excitation of the He I Lines 


±0 


0005 


Reddening correction 


±0 


0005 


Optical Depth of the He I Triplet Lines 


±0 


0003 


He I and H I Line Intensities 


±0 


0003 



4. The Four Main Sources of Error 

The most important source of error is the collisional excitation of Balmer lines. 
Neutral hydrogen atoms in excited states may form not only by the usual pro- 
cess of H + recombination, but also through collisions of neutral hydrogen with 
electrons. The recombination cascade that follows such excitations contributes 
to the observed intensity of Balmer lines, mimicking a larger relative hydrogen 
abundance and, hence, a smaller helium abundance. To estimate this contri- 
bution it is necessary to have a tailored photoionization model for each object 
that fits properly the temperature structure. The contribution to the Balmer 
line intensities depends strongly on the temperature: therefore this effect, and 
consequently the associated error in its estimate, increases for H n regions of 
lower metallicity and consequently higher temperature. 

The second most important source of error is the temperature structure. 
Most determinations neglect the possible presence of temperature variations 
across the H n region structure and assume that T(0 in) is representative of 
the zone where the He I recombination lines form. However, other temperature 
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determinations based on different diagnostics yield lower values; furthermore, 
photoionization models do not predict the high T(0 in) values observed. These 
results indicate that temperature variations are indeed p resent in H II regions , 
and this result should be included in the Y determination ( Peimbert et al. 2002). 
The best procedure to take into account the temperature structure is to self- 
consistently determine T(He n) from a set of He I lines by means of the maximum 
likelihood method. The Y abundances derived from T(He n) are typically lower 
by about 0.0030 — 0.0050 than those derived from T(0 in). The difference 
between both temperatures does not have a significant trend with the metallicity 
of the H II region, hence the systematic error introduced by the use of T(0 in) in 
the Y determination is similar for objects with different metallicities. The error 
quoted under "Temperature Structure" in Table 1 is the resid ual error due to 
the un certainty in the T(He Ii) determinations of the dataset bv lPeimbert et aT 
(120071 1. 

The third most important source of error is the extrapolation of the derived 
Y values to zero metallicity through the AY/ AO ratio. This problem will be 
discussed in the next section. The fourth most important source of error is the 
uncertainty on the recombination coefficients of the He I lines. 



5. AY /AO from Models and Observations 



To determine the Y p value it is customary to use the Y values of a set of O poor 
galaxies and to extrapolate the Y values to the case of O = using the following 
equation: 



Y„ = Y-0 



AY 

AO' 



(1) 



where O is the oxygen abundance per unit mass. To obtain an accurate Y p value, 
a reliable determination of AY /A O for O-poor ob j ects i s needed. 

The AY/ AO value derived bv lPeimbert et ail ( 2000l ) from observational re- 
sults and models of chemica l evolution of ga laxies amounts to 3.5 ± 0.9. More 
recent results are those by iPeimbertl (12003] ) who finds 2 .93 ± 0.85 from ob- 
servations of 30 Dor and NGC 346, and by llzotov et al. ( 20061 ) who, from the 
observations of 82 H n regions, find AY/ AO = 4.3 ± 0.7. We have recomputed 
the value by Izotov et al. considering two systematic effects not considered by 
them: the fraction of oxygen trapped in dust grains, which we estimate to be 
10% for objects of low metallicity, and the increase in the O abundances due 
to explicit taking into account the presence of temperature fluctu ations, which 



for t his type of H n regions we estimate to be about 0.08 dex (jRelaho et al 
2002). From these considerations we obtain for the Izotov et al. sample a 
AY /AO = 3.2 ± 0.6. 

From chemical evolution models of galaxies it is found that AY/ AO depends 
on: the stellar yields, the initial mass function, the star formation rate, the age, 
and the O value of the galaxy in question. Models with substantial outflows of 
O-rich material can produce large AY/ AO rat ios but they are ruled out by t he 
low C/O values observed in irreg ular galaxies (ICarid et al.lll995l . Il999l . l200d ). 

Carigi &: Peimbert ( 20071 ) have produced chemical evolution models of the 
following types: closed box, inflow of gas, and outflow of gas of well-mixed 
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Figure 1. Two chemical evolution models for NGC 6822 with different gas 
infall and outflow histories, and the same star formation rate derived from 
observations. The mo dels are the 7L and 8S (continuous and dashed lines, re- 
spectively) studied by ICarigi et al.l (|2006f ) . The first panel shows the increase 
of the helium and oxygen abundances relative to the primordial values, the 
second panel shows the gaseous content as a function of time, which is widely 
different for the two models, and the third panel shows the star formation 
history, which is the same one in both models. 



material. They find that AY/ AO is practically constant for models with the 
same IMF, the same age, the same star formation history, and an O abundance 
smaller than ~ 4xl0 -3 . They find also that 2.4 < AY/ AO < 4.1 for models 
with different star formation histories and different values of the upper mass 
limit of the IMF. The results derived from the chemical evolution models are in 

very good agreement with the observations mentioned a bove. 

Based on the observations and models discussed, iPeimbert et al.l ( 20071 ) 



adopted a value of AY/ AO = 3.3 ± 0.7 in the computation of the error budget 
presented in Table 1 and the Y p value presented in Table 2. 
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For models of galactic chemical evolution that reach values of O > 4 x 10~ 3 
at present time, the AY/ AO ratio of the interstellar medium increases with the 
O abundance due to two effects: the helium production by low and intermediate 
mass stars and the increase of the helium yield of massive stars due to stellar 
winds. 

As an example of the minor role that well mixed outflows play in the 
chemical evolution of galaxies in Figure 1 we present the AY versus AO be- 
havior for two chemical evolution models of NGC 6822 ( Carigi et al.l 120061 ; 
Carigi Sz Peimbert 20071 ). These models have the same SFH, which was de- 
rived from observations, but are drastically different in their gas flow histories 
and show practically the same AY/ AO behavior. 

To compute stellar evolution models with O < 4 x 10~ 3 we propose to use 



the following relation for the initial Y and O abundances 
Y = 0.2474 ± 0.0028 + (3.3 ± 0.7)O. 



(2) 



For O > 4 x 10~ 3 (ICarigi k Peimberd[2007h . Y increases faster with the increase 
of O than in the previous equation. 

The ratio AY/ AO can also be expressed in terms of AY/ AZ if the fraction 
of O/Z is known. Based on observations of galactic and extragalactic H n 
regions, we propose to assume that, for models with O < 4 x 10~ 3 , O constitutes 
55% ± 5% of the total Z value, implying AY / AZ = 2.0 ±0.6. For O > 4 x 10" 3 
the fraction of Z due to O decreases due to the increase of the C/O, N/O 
and Fe/O ratios with the increase of O. From the chemical composition of the 
Orion nebula and M17 it is found that t he fraction of Z due to O drops to 
about 42% for 0.0057 < O < 0.0082 (e.g. iPeimbertl 120031 : lEsteban et all 12004 : 
Garcfa-Roias et al.ll2007l . and references therein). 

Previous determinations of the AY I AZ ratio based on models or observa- 
tions have been in the 1.0 to 6.0 range (e.g. |Peimberd[~L995l ; Fukugita &: Kawasaki 

2006, and references therein). It is interesting to note that 25 years ago lChiosi &: Matteucci 
(| 19821 ) determi ned a value of AY / AZ ~ 2.0, in excellent agreement with the re- 
cent results by ICarigi fc Peimbertl ||2007l ). 



6. Discussion 



In Table 2 we present some of the best Y p determinations of the last few years. 
The Yp values and their statistical errors are the ones presented in the original 
papers. After the statistical error we list the systematic error estimated by us, 
which depends on one or more of the following sources: a) the change in the 
published emissivities of the He I lines ( Porter et al. 20051 ) ; b) the change in the 
publi shed collisional excitation coefficients of the H I lines ( Anderson et al. 20001 . 
2002); and c) the temperature structure of the H II region. Adopting the new 
He I emissivities, the Y determination of individual H n regions is increased by 
about 0.004. T he change in the H I c ollisional excitation coefficients goes in the 



same direction ( Peimbert et al. 20071 ). although in this case the size of this effect 



varies from object to object and a tailored photoionization model for each object 
is needed to obtain a good estimate of it. Both effects pr oduce an increase in th e 
Y p determination: for the sample of H n regions used bv lPeimbert et all ( 20071 ). 
the increase in Y p due to the adoption of the new He I emissivities amounts to 
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about 0.0040, while the increase due to the adoption of the new H I collisional 
excitation coefficients amounts to 0.0025. Finally, the Y value of individual H n 
regions decreases when the temperature structure of the H II region is taken 
into account, since in the self-consistent solutions for all the observed He I line 
intensities, the lower T(He n) values imply higher densities, the higher densities 
produce a higher collisional contribution to the He I intensities, a nd consequently 
l ower heliu m abundances. For th e objects in the sample used bv lPeimbert et all 
(2002) and Peimbert et al. ( 20071 ). this effect decreases the Y determinations by 

amounts ranging from 0.003 to 0.009. 

The Y p determination by llzotov et al.l ( 19991) is affected by all thr ee of 
the above sourc e s of systematic error; those by iLuridiana et al.l (2003) and 
Izotov fc Thuan J 2004! ) are affected by sources a) and b), while the one by 



Izotov et al.l (120061') is affected by so urce b). The disagreement between the 



Y p derived by Luridiana et al.l ( 20031 ) with the Y p derived from WMAP under 
the assumption of SBBN implies the need for "new physics" . The new physics 
ne eded to reconcile the two Y v v alue s turned out to be th e "new atomic physics" 
bv I Anderson et all (120001 . I2002T) and lPorter et al.1 (j2005T ). With the new physics 
Peimbert et al.1 (|2007l ) found agreement, within the observational errors, be- 



tween the observed Y p value and that derived from the WMAP results under 
the assumption of SBBN. 

The Y p derived bv lPeimbert et all (|2007l ) together with SBBN implies that 
tt b h 2 = 0.02054 ± 0.00639 (jSteigman! l2006al lbh . where tt b is the baryon closure 
parameter, and h is the H ubble paramet e r. Th is value is in excellent agreement 
with the value derived bv ISpergel et al. ( 20061 ) from the WMAP results under 
the assumption of SBBN, which amounts t o Q. h h 2 = 0.0223 3 ± 0.0 0082. 

The c omparison of the Y p derived by iPeimbert et al. (2007) with the Y p 
derived by ISpergel et al.l ( 20061 ) from the WMAP data together with the as- 
sumption of SBBN provides strong constraints for the study of non SBBN (e.g. 
Cvburt et al.ll2005l ; ICoc et al.ll2006l . and references therein). 



In Table 2 we also present our Y p prediction for 2010. We consider that in 
the next few years it will be possible to reduce the statistical errors in the Y p 
determination to about 0.0020 by obtaining a new set of observations of brighter 
and slightly O-richer H n regions than the ones that have been used so far. A 
more extensive discussion of the relative advant ages of these H I I regio ns with 
respect to more metal-poor ones can be found in Peimbert et al. ( 20071 ) . 



7. Summary and conclusions 

In this contribution we have presented some recent determinations of Y p and 
discussed the reasons underlying the differences betwe en them. For the most 
recent determination, the one bv lPeimbert et ah ( 20071 ). we have presented the 



error budget in terms of thirteen different sources of error. The AY/ AO ratio 
which enters as a crucial factor in one of the three main sources of error, has 
been discussed in the light of recent observations and models of galactic chemical 
evolution. 

The Y p determinations by different groups are slowly converging among 
them as systematic errors are progressively identified and correct ed for. The need 



for "new physics" that had been suggested by recent results (e.g. . ILuridiana et al 
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Table 2. Primordial helium abundance values a 



Source Y p 

Izotov et al. (1999), this work 0.2452 ± 0.0015 ± 0.0100 

Luridiana et al. (2003), this work 0.2391 ± 0.0020 ± 0.0070 

Izotov & Thuan (2004), this work 0.2421 ± 0.0021 ± 0.0075 

Izotov et el. (2006), this work 0.2462 ± 0.0025 ± 0.0040 

Peimbert et al. (2007) 0.2474 ± 0.0028 

Prediction (2010), this work 0.2??? ± 0.0020 

Spergel et al. (2006) 0.2482 ± 0.0004 

a Direct Y p determinations based on observations of H n regions, with the exception 
of that by Spergel et al. (2006), which is based on the baryon to photon ratio derived 
from WMAP and the assumption of SBBN. 



2003) seems now to b e fulfil led b y new atomic physi cs, i.e. the atomic data by 
Anderson etaD (|2000l . I2002T ) and iPorter et al.1 (120051 ). On the other hand, the 



temperature structure of H n regions is still a source of systematic error if an 
appropriate scheme for temperature is not adopted. The proper temperature to 
determine the helium abundance is T(He n), derived self-consistently from the 
intensities of the helium lines. Adopting this temperature, the Y p value derived 
from H II regions agrees with the Y p derived from the WMAP data under the 
assumption of SBBN. On the other hand, the use of T(0 in) to determine the Y 
values from H n regions produces Y p values from 0.003 to 0.006 higher than those 
found adopting T(He n), that is Y p values more than la higher than the one 
predicted by the WMAP observations combined with the assumption of SBBN. 

It is important to continue the effort on the study of the physical conditions 
in H II regions, this effort will permit us to lower the error on the Y p determi- 
nation, which in turn will permit us to improve our knowledge on the possible 
importance of non SBBN. 
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